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ABSTRACT 


Rate  equations  and  enthalpic  parameters  were  derived  for  use  in  differential  thermal 
analysis  (DTA)  of  polymerization  reactions*  General  free  radical  polymerization  kinetics 
are  discussed  and  the  existence  of  a  "stationary  state"  under  DTA  conditions  was  demon¬ 
strated,  by  comparison  of  a  theoretical  computer  curve,  based  on  stationary  state  assump¬ 
tion,  with  actual  experimental  conditions.  These  curves  appear  to  be  almost  superimposable. 

This  technical  documentary  report  has  been  reviewed  and  is  approved* 
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INTRODUCTION 


Differential  thermal  analysis  techniques  have  in  recent  years  been  extended  to  encompass 
the  field  of  solution  reactions.  Borchardt  and  Daniels  (Reference  1)  were  successful  in 
measuring  the  heat  of  reaction  and  in  determining  the  kinetic  parameters  for  the  decomposi¬ 
tion  of  benzene  diazonium  chloride  and  the  reaction  of  N,  N-dimethyl  aniline  and  ethyl  iodide 
in  solution.  This  technique  has  now  been  adapted  for  use  in  determining  polymerization 
kinetics  and  thermodynamic  parameters.  In  Part  I  of  this  report  (Reference  2)  the  apparatus 
and  necessary  modifications  are  described.  In  the  present  report  the  derivation  of  the  rate 
expression  and  enthalpy  relationship  are  presented,  which  permit  the  extension  of  Borchardt 
and  Daniels*  work  to  polymerization  systems. 

Borchardt  and  Daniels  (Reference  1)  have  given  the  scheme  by  which  rate  of  reaction  and 
enthalpy  data  can  be  obtained  by  differential  thermal  methods  from  simple  first-order  and 
pseudo-first-order  reactions  taking  place  in  solution.  Starting  with  the  assumption  that  in 
a  dilute  solution  the  heat  capacities  (at  constant  pressure)  of  solvent  and  that  of  solvent  plus 
reactants  are  very  nearly  the  same,  the  heat  generated  or  absorbed  by  a  chemical  reaction 
is  given  by: 

dH  =  Cp  dAr  +  kAT  dt  (1) 

where  K  is  the  cell  constant,  and  A  T  is  the  height  of  the  curve  in  degrees  C.  Integration 
of  this  equation  between  t  =  0  and  t  =  oo  gives  the  total  heat  transferred,  A  H, 

flO 

Ah  =  Cp  (  AToo  -  ATq)  +  kJ  At  dt  (2) 

0 

Since  AT  is  zero  both  at  t  =  0  and  t  =  00  ,  the  first  term  is  also  zero,  reducing  the  equation 
to: 


Ah  =  KA  (3) 

A  being  the  total  area  under  the  curve.  This  condition  expresses  then  the  enthalpy  of  the 
reaction.  A  typical  DTA  curve  is  shown  in  Figure  1. 


Figure  1.  Typical  DTA  Curve 


Manuscript  released  by  the  authors  21  April  1964  for  publication  as  an  RTD  Technical 
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Assiuning  that  the  heat  evolved  during  a  small  time  interval  is  directly  proportional  to 
the  number  of  moles  reacting  durir^  that  time: 

dH  =  dn  ,  (4) 

=  total  number  of  moles  that  react),  then  substitution  of  equation  (4)  into  equation  (1) 

gives: 


M  ^  11  (  r  y 

dt  KA  I  dt  *  ) 


dAT 

where  -g^ —  is  the  slope  and  AT  the  height  of  the  curve  (Figure  1).  The  number  of  moles 
unreacted  (n)  at  any  instant  is  equal  to  the  initial  moles  (n^)  minus  the  number  of  molejs 
that  have  reacted: 


0 


Substituting  equation  (5)  into  equation  (6)  and  integrating  gives: 


n  =  no  -  C  CpAl  +  Ko  )  (7) 

The  quantity  a  is  shown  in  Figure  1.  The  e^q^ression  for  the  rate  constant  of  a  reaction  of 
order  x  with  respect  to  one  component  is: 


dn  ^  kn 
dt  ^ 


ond  k  =  — 


V  dt 


where  k  =  reaction  rate  constant,  V  =  volume,  and  n  =  moles.  If  n^  =  then  substitution 
of  equations  (5)  and  (7)  into  equation  (8),  and  rearranging  gives: 


'  KAV  1 

.  J 


Cp  +  kAt 

K(A-a)  -  Cp  At] 


For  a  first  order  reaction  this  becomes: 


K(A-a)  -  CpAr 


The  relative  magnitude  of  the  terms  in  equation  (10)  is  such  that  those  involving  are 

usually  much  smaller  than  the  quantitites  from  which  they  are  added  or  substracted.  By 
neglecting  Cp*  equation  (10)  may,  therefore,  be  simplified  to: 


2 


k 


At 

A-o 


(11) 


This  is  essentially  the  equation  used  by  Borchardt  and  Daniels,  and  their  kinetic  results  are 
in  good  agreement  with  data  obtained  by  conventional  methods.  They  also  obtained  enthalpy 
data  on  the  reactions  mentioned,  but  no  comparisons  are  possible  since  no  literature  values 
are  available. 

These  same  authors  further  extended  their  derivations  to  include  reactions  between  several 
components  of  the  form  iL  +  mM  - ►  products,  for  which  the  rate  expression  becomes: 

=  (L-z/  (M  -  z)*"  (12) 


in  which  case  equation  (10)  is  expanded  to 


_  +  m-1  _  . 

[^] 

[kCa-o)  -  CpAr]  [(K^  A  - 0 ) -CpAr]' 


POLYMERIZATION  KINETICS 

The  above  formalization  is  based  upon  stoichiometric  reactions  and  cannot  as  such  be 
used  to  describe  free  radical  polymerization.  It  will  be  necessary,  therefore,  to  develop 
the  required  relationships  for  these  systems.  The  series  of  reactions  which  make  up  the 
overall  process  of  free  radical  polymerization  may  be  written  as: 

1.  Generation  of  Initiator  Fragments:  I^ — ^  21*; 

2.  Chain  Initiation:  I*  +  M  — ►  M-  ;  [I.]  [M] 

3.  Chain  Propagation:  +  M  — ►  M  -j+j  ;  Rp  =Kp[M][M] 

4  .  Chain  Termination:  IVI*j  +  M-j  - ►  Inactive  Products:  =  2k^ 

where  f  is  the  efficiency  of  initiator  fragments.  Because  of  lack  of  knowledge  of  the  con¬ 
ditions  in  this  experiment,  f  was  assumed  to  be  1. 

The  overall  rate  of  polymerization,  the  disappearance  of  monomer,  or  the  formation  of 
polymer  is  given  by: 


If  the  kinetic  chain  length  is  reasonably  high,  the  loss  of  monomer  by  the  chain  initiation 
step  is  negligible  compared  with  that  lost  during  propagation  and  the  rate  of  polymerization 
may  be  taken  as  equal  to  the  rate  of  propagation: 


-  d  [m] 
dt 


kp  [M-]  [m] 
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In  general,  the  heat  evolved  or  absorbed  during  polymerization  is  made  up  of  the  heats 
liberated  or  absorbed  by  each  of  the  several  reactions  which  take  place: 

^  (  A  +  A  H|  +  A  Hp  +  A  +  A  Hgt  + . )  , 


where  A  Hfji  is  the  total  heat  and  the  subscripts  refer  to  initiator  decomposition,  chain 

initiation,  propagation,  termination,  and  chain  transfer,  respectively.  Even  though  the 
heat  generated  or  absorbed  per  mole  for  all  these  reactions  may  be  significant,  due  to 
the  fact  that  the  number  of  initiator  decomposition,  chain  initiation,  chain  termination, 
and  chain  transfer  reactions  that  occur  per  propagation  step  is  very  small  (about  equal 
to  the  reciprocal  of  the  number  average  degree  of  polymerization)  the  heat  generated 
during  polymerization  is  with  very  small  error  due  to  the  propagation  reaction  alone. 

Equation  (14)  contains  two  easily  measured  quantities^  -d  [M  ]/dt  and  [M],  one  quantity 
which  is  difficult  to  measure,  LM‘],  and  one  quantity,  the  rate  constant,  |the  determina¬ 
tion  of  which  is  one  of  the  objectives  of  this  work.  Free  radical  concentrations  generally 
are  quite  small,  in  the  10^^  to  10-^^  region,  and  no  accurate  methods  are  available  for 
their  measurement.  It  is  possible  to  circumvent  the  measurement  of  this  quantity  by  use 
of  the  stationary  state  assumption.  This  assumption  may  be  stated  as:  the  rate  of  change 
of  radical  concentration  is  very  much  smaller  than  either  its  rate  of  formation  or  rate  of 
disappearance.  Since  the  rate  of  formation  of[M]  is: 


+  d  [m  •] 

dt 


(15) 


and  its  rate  of  disappearance  is: 

-  d  [m  •] 
dt 


2kt  [m.]' 


(16) 


the  instantaneous  M*  concentration  is: 


d  [m-] 

dt 


k,  [I.]  [M]-2kj[M.]'' 


(17) 


According  to  the  stationary  state  assumption  d  [M*  ]/dt  is  very  much  smaller  than  either 
kj  [l-]  [M]or2kt[M-]^,  therefore: 

k,  [!•]  [m]  S  2k^  [m.]^ 

and 

[m-]^  kj  [i  ]  [m]  (18) 

A  similar  stationary  state  assumption  may  be  used  to  replace  the  I*  in  the  expression 
d[l*]  /dt  =  2fkd  [iJ-kj  [!•]  [m]  itherefore 

2fkj  [ig]  =  [I*]  M  (19) 
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Combination  of  equations  (18)  and  19)  and  substitution  into  equation  (14)  lead  to: 


d  [m] 

"dt 


Pa] 

kt 


(20) 


The  validity  of  the  stationary  state  assumption  has  been  borne  out  by  numerous  experi¬ 
ments  conducted  at  constant  temperature.  It  has  been  subjected  recently  to  a  critical  review 
and  the  general  criterion  has  been  derived  for  its  use  in  isothermal  experiments 
(Reference  3).  However,  little  is  known  of  the  validity  of  the  assumption  under  differential 
thermal  analysis  conditions  where  the  temperature  changes  continually  throughout  the  re¬ 
action.  It  is, therefore,  necessary  to  examine  this  in  some  detail. 

CONSIDERATION  OF  “STATIONARY  STATE  ASSUMPTION”  UNDER 
CONSTANTLY  INCREASING  TEMPERATURE  CONDITIONS 

For  the  investigation  of  the  applicability  of  “stationary  state”  under  conditions  here, 
the  polymerization  of  styrene  initiated  by  AIBN  was  selected,  since  its  rate  constants 
are  better  known  and  the  initial  and  final  concentrations  and  the  heating  rate  could  be 
experimentally  determined.  Values  for  kp  and  k^  were  taken  from  the  literature  (Refer¬ 
ence  4).  Initiator  decomposition  parameters  were  determined  at  this  laboratory.  Five 
values  covering  the  range  of  10®  to  10®  were  assigned  to  Iq,  assuming  it  to  be  of  the  same 

general  order  of  magnitude  as  the  propagation  rate  constants.  These  data  are  summarized 
in  Table  1. 


TABLE  I 


Rate  Constants  in  Styrene  Polymerization 


Reaction 

Rate  Constant 

A. 

1 

E.,  cal. 

1 

Initiator  decomposition 

6.2  X  10^® 

32,000 

Termination 

*^t 

1.32  X  10® 

2,400 

Propagation 

k 

P 

2.2  X  10’ 

7,800 

initiation 

Vl 

2.2  X  10’ 

7,800 

■^1-2 

CD 

O 

10,000 

CO 

1 

(0 

O 

«-4 

6,000 

k.  . 

1-4 

10® 

10,000 

k.  - 
1-5 

10® 

20,000 

The  initial  initiator  concentration  I^  was  5.25  x  10“®  mole/1,  the  initial  monomer  con¬ 
centration  M^  =  1.5  mole/1,  the  heating  rate  dT/dt  =  1.027  x  10"®  *’K/sec.  The  reaction 
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visibly  started  at  =  341.5'’K  and  ended  ('^  )  at  390“K  (approx,  duration  4800  sec); 

T_  -  T  =  T„  is  dierefore  48.5*K  and 
00  o  R 


TABLE  2 


This  approximation  results  in  a  considerable  simplification  of  subsequent  calculations, 
thereby  more  than  offsetting  the  small  error  introduced. 

Accordingly, 

-It}”  rT^}  (23) 
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(24) 


Since  =  (dT/dt)t,  the  rate  constants  can  be  expressed  in  the  following  form: 
K 


.  «  h:t 

k.  =  Bje  • 


^  where 


B.  = 


Aje 


RTn 


,  ond  hj  = 


dt 

RTn* 


and  h^  are  constants  for  a  particular  experiment  and  their  values  may  vary  depending 
on  the  heating  rate,  dT/dt,  and  the  initial  temperature  of  the  reaction,  T^, 

The  rate  constants  were  then  calculated  and  are  given  in  Table  3. 


TABLE  3 


Rate  Constants  -  Approximate  Form 


Rate  Constants 

B. 

1 

h.  X  10^ 

1 

2.8  X  10-® 

14.0 

3.9  X  10'' 

1.06 

k 

P 

2.4  X  10® 

3.43 

^i-1 

2.4  X  10® 

3.43 

’^i-2 

4.8  X  10-^ 

4.4 

^i-3 

1.5  X  10® 

2.6 

k. 

1-4 

4.8  X  10^ 

4.4 

k.  . 

1-5 

2.9  X  10-® 

1.76 

Since  - 


] 


dt 


=  [I3  ],  integration  of  this  equation  gives: 

[I2]  =  lo  exp{-/  k(jdt]>  •, 

0 


denoting  2  fk^j  [I2  ]  as  g,  we  obtaih: 


g  -  Stk^IpC 


-/  V 
0 


(25) 
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The  exact  rate  expressions  were  calculated: 


=  3.0  X 


f  e%p 


4  X  l< 


t  -  2  X 


1.4  n  l( 


(26) 


dt 


=  g-kj 


(27) 


d  [m*] 
dt 


7.8  X 


X  e 


1.06  X  10  t 


(28) 


d  [m] 
dt 


-2.4  X 


X  e 


3.43  X  10  t 


(29) 


For  the  approximate  equations  using  the  stationary  state  assumption  the  following 
equation  was  used  in  lieu  of  equations  (27)  and  (28): 

g  =  7.8  X  lo"  X  ^  [m-]'  =  k,  [i.][m]  (30) 

These  two  sets  of  equations,  (a)  (26-29)  and  (b)  (26),  (30)  and  (29)  were  solved  independ¬ 
ently  for  k.  by  digital  computer.  For  purposes  of  cilchlation  f  was  assumed  to  be  1, 

’  ^  ,  d  [m] 

Figure  2  shows  the  functions  g,  [M*3,  ,  [l*  ],  (fork,  j),  and[rg] 

(for  kj  g)  as  calculated  by  computer.  The  fimction  d  CM]/dt  was  calculated  from 

data  obtained  by  computer  for  the  function  [M].  Figure  3  describes  the  functions  [M] 
and  [M*]  showing  both  exact  as  well  as  approximate  equations. 

The  functions  obtained  from  the  two  sets  of  equations  were  practically  identical.  The 
maximum  deviation  of  the  function  [M]  was  approximately  0.15  percent  and  that  of  [M*], 
0.4  percent,  the  mean  deviations  being  negligible.  The  deviations  of  [T  ]  calculated  from 
both  sets  of  equations  and  their  effect  on  [M*]  and  [M]  were  so  small  fliat  they  could  not 
be  detected  by  the  8-dlglt  computer  and  may,  therefore,  also  be  considered  negligible. 

The  plot  of  the  functions  [M*  ]  and  -d  [M]/dt  shows  an  abrupt  Increase  shortly  after  the 
start  of  the  reaction  coinciding  thereafter  closely  with  the  curve  derived  from  the  station¬ 
ary  state  assumption.  The  evaluation  of  the  rate  constants  from  experimental  curves 
will  have  to  take  this  into  account  by  curve  fitting  in  this  region  and  adjustment  of  the 
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base  line.  This  condition  is  shown  in  Figure  4  by  the  slanted  base  line  of  the  experimental 
curve. 


Assuming  proportionality  of  AT  and  [M]  as  well  as  the  adjusted  base  line,  the  experi¬ 
mental  curves  were  plotted  by  means  of  the  least  squares  method.  The  resulting  curve 
was  superimposed  on  the  theoretical  curve  of  [M]  (Figure  4).  Except  for  the  afore¬ 
mentioned  abrupt  shift  at  the  beginning  of  the  reaction  there  is  good  agreement  between 
the  two  curves.  The  deviation  of  the  experimental  curve  from  the  calculated  values  in¬ 
dicates  that  the  reaction  must  have  started  in  a  very  gradual  manner  before  it  reached 
the  point  indicated  by  the  thermocouple  as  the  start  of  the  reaction.  This  analysis  has, 
nevertheless,  shown  that  the  stationary  state  assumption  is  quite  valid  within  the  general 
conditions  of  the  experiment. 

DERIVATION  OF  RATE  EQUATION 


By  neglecting  the  C 


P 


te  rm,  as  before,  equation  (5)  may  be  reduced  to: 


d  n  . 

dt  ■  A 


(31) 


The  decomposition  of  the  initiator  which  essentially  goes  to  completion  (i.e.,n^  =  I^) 
may  be  stated  as: 


flpAT' 

VA' 


(32) 


where  T’  and  A’  are  used  to  distinguish  the  initiator  contribution  from  the  rest  of  the  re¬ 
action  terms.  For  practical  purposes  T’  and  A’  are  determined  in  a  separate  run  of 
initiator  decomposition  conducted  under  the  same  conditions  as  the  polymerization  reaction. 

Since  the  disappearance  of  monomer  is  equal  to  the  formation  of  pol5nner: 

_  d  M  .  [P]  .  pAi 

dt  ■  dt  ■  VA 

where  P  is  the  total  amount  of  polymer  formed.  Similarly  to  equations  (6)  and  (7),  the 
monomer  concentration  at  any  given  point  is  again  equal  to  its  initial  concentration,  M^, 

less  the  amount  of  polymer  formed: 


M  =  (34) 

It  has  already  been  pointed  out  (equation  14),  that  the  propagation  rate,  now  expressed  as 
dP/dt  may  be  assumed  to  be  equal  to  the  overall  rate, 

.  K,  [M]  [M.] 
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Rearrangement  of  equation  (20)  and  substitution  of  equations  (33)  and  (34)  then  gives: 


I 


_ pAt _ 

(MqA-  Pa)  Hkj  j^lgj 


I 


(35) 


Substitution  of  equation  (32)  into  the  above  equation  gives  the  final  expression  for  the 
polymerization  rate  constant: 


^  _ cAt _ 

(MqA-po)  (  At' 

where  the  constant  C  is: 


(36) 


I 


Reintroduction  of  the  heat  capacity  term  expands  equation  (36)  to: 


(37) 


kp  _  CiCp  ^  +  kAt) 

k^^  (  At' )^  [k(MoA-po)  -  PCpArJ 


(38) 


UNITS  AND  DIMENSIONS 

The  reaction  rate  is  in  liters,  mole  ,  sec  ;  M^,  I^,  and  P  in  moles;  V  in 

liters;  A,  A’  and  a  in  “Csec,  AT  and  AT’  in  °C;  in  calories  dAT/dt  in  '’Csec~^; 

and  K  in  calories  sec . 

The  cell  constant  K  needs  only  be  evaluated  once  from  a  reaction  with  known  enthalpy. 
C  is  the  total  heat  capacity  of  the  solvent  and  monomer,  but  since  we  are  dealing  with  a 

Jr 

dilute  system  the  heat  capacity  of  the  solvent  is  usually  assumed  to  be  equal  to  that  of  the 
system,  d  AT/dt  is  the  slope  of  the  DTA  curve.  The  enthalpy  is  in  calories  mole”’^ . 

The  rate  constants  are  plotted  in  the  usual  manner  as  -log  k  vs.  l/T^K  and  from  the 

activation  energy  plot  the  activation  energy  E  -  E  /2  and  the  frequency  factor  A  /A 

p  t  p  t 

can  then  be  obtained. 

This  work  demonstrates  the  feasibility  of  extending  DTA  into  the  realm  of  polymeriza¬ 
tion  kinetics.  It  has  been  shown  that  stationary  state  conditions  prevail  when  reaction  is 
subjected  to  a  moderate  heating  rate.  The  present  heating  capacity  of  the  apparatus  lies 
between  approximately  0.  5  and  2.5’’/mln  and  while  it  is  to  be  expected  that  some  increase 
in  the  heating  rate  beyond  this  limit  should  upset  the  stationary  state  the  exact  practical 
limits  are  still  to  be  determined. 
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The  relative  ease  with  which  kinetic  parameters  as  well  as  enthalpy  can  be  determined 
by  means  of  DTA  reduces  the  amount  of  work  involved  from  that  extending  over  a  period 
of  weeks  and  even  months  to  one  day. 
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SCALE  OF  FUNCTIONS 
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from  co; 


OEXPERIMENTAL  DTA  CURVE 


Figure  4,  The  theoretical  and  experimental  function 


